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FATIGUE LIFE ANALYSIS FOR TRACTION DRIVES WITH 
APPLICATION TO A TOROIDAL TYPE GEOMETRY 
by John J. Coy, Stuart H. Loewenthal, and Erwin V. Zaretsky 

Lewis Research Center and 
U. S. Army Air Mobility R&D Laboratory 

SUMMARY 

A contact fatigue life analysis for traction drives has been developed based on a 
modified Lundberg-Palmgren fatigue theory which is the accepted basis of rating rolling- 
element bearing life in the bearing industry. The analysis was used to predict the ser- 
vice life of a modern cone-roller toroidal traction drive. Material and design life ad- 
justment factors were considered. A constant operational traction coefficient of 0.06 
was chosen. A minimum contact overload factor of 1.2 was selected to preclude roller 
slip. An elastohydrodynamic film thickness calculation was performed for the input and 
output contacts at a representative duty cycle operating point. Life adjustment factors 
due to advances in materials and lubricant technology were considered. 

The traction drive geometry consisted of a double roller per cavity dual cavity 
toroidal drive with a torus diameter of 11.43 centimeters (4.50 in.), a cavity diameter 
of 8.84 centimeters (3.48 in.) and a roller crown radius of 3.12 centimeters (1.23 in.). 
Parametric results were presented for life and Hertz contact stress as functions of load, 
drive ratio, and size. 

To determine vehicle service life, the toroidal drive was incorporated into a typical 
compact passenger vehicle drive train using a representative duty cycle condition of 
29 engine horsepower, 48.3 kilometers per hour (30 mph) vehicle speed, and a drive 
ratio of 1.2. The mileage life was calculated at 19 200 kilometers (11 900 miles). This 
is 38 times less than previously stated in the literature for the same condition. At this 
nominal operatii^ condition the minimum lubricant film thickness to composite surface 
roughness ratio was calcinated as 0. 56 on the output contact. The film thickness is 
marginal to prevent severe roller surface distress. 

At severe operating conditions (high drive reduction ratios and high transmitted 
loads) , the brinelling limit for the input disk -roller contact is approached . Shock or 
dynamic overload protection is marginal. 

The present fatigue analysis has shown that life is proportional to the 8 . 4 power of 
drive size. To provide for 160 900 kilometers (100 000 miles) of vehicle service at the 
aforementioned duty cycle point, the cone roller toroidal drive size must be increased to 
1. 3 times its current size. 



INTRODUCTION 


With the increased concern for the nation's dwindling energy resources, the pas- 
senger automobile has become a prime target for energy conservation efforts . As an 
example, in 1970 the transportation sector accounted for more than 24 percent of the 
total energy consumed in the U.S. and the gasoline consumption of private motor ve- 
hicles represents about 55 percent of this portion (ref. 1) . 

Of the several areas where improvements in vehicle fuel economy can be realized, 
altering the transmission system is one of the most rewarding (ref. 2). Several in- 
dependent studies have estimated that the incorporation of a continuously variable trans- 
mission (CVT) in place of the conventional automatic transmission in a standard piston 
engine vehicle would result in a fuel mileage gain of up to 28 percent over the federal 
driving cycle (refs. 2 and 3). Furthermore, the use of a CVT in conjunction with other 
vehicle and engine improvements for future vehicles could mean as much as a 60 per- 
cent fuel savings (ref. 2). These savings are attributed to two benefits of the CVT. 

With proper controls the CVT is able to regulate the ermine speed to operate along a 
minimum fuel consmnption line. Thus, fuel economy is improved. Also, since better 
vehicle acceleration is obtained with the CVT through engine speed control, smaller 
engines may be used adding to the fuel economy . 

The goal of developing a commercially viable automotive CVT has brov^ht about 
renewed interest in traction type transmissions. These transmissions transmit power 
through shear of the lubricant film at the surface of smooth roller elements which are 
in highly loaded contact. Speed ratio changes are accomplished by the repositioning of 
some of the drive rollers which in turn alters the contact rolling radii of the components 
in the system. The principal advantages that traction CVT’s hold over hydraulic types 
is in their ability to transmit power smoothly and quietly with high efficiency, normally 
greater than 90 percent. 

The automotive traction transmission made its first commercial appearance in the 
1909 Carter car. In the 1930's General Motors introduced the Transitorq (ref. 4) and 
Austin produced the Hayes transmission (ref. 5). Neither of these toroidal traction 
transmissions were in service for very long. One factor that hindered the widespread 
acceptance of the earlier traction drives was the limited durability of roller component 
materials and the imcertainties in life prediction methods. High quality roller materials 
and lubricants are essential for long service life since high contact pressures often 
greater than 240 000 newtons per square centimeter (350 000 psi) maximum Hertz stress 
are needed to transmit significant power throi^h these drives. Since the drive's design 
life is limited by roller fatigue which is a strong function of contact pressure, high 
powered traction drives have normally been restricted to limited-life applications 
(ref. 4). However, the power ratings and service life of modern traction drives have 
improved somewhat due to the production of cleaner bearing steels with better fatigue 
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resistance and the development of lubricants which possess better tractive properties 
(ref, 6), 

Recently, an advanced toroidal traction transmission has been proposed for auto- 
motive application (ref. 7) . Use of this CVT could greatly improve vehicle fuel econ- 
omy but its acceptance depends heavily on conversion cost factors and its reliability for 
long term vehicle service. 

Present methods to predict the design lives of traction elements are limited to em- 
pirical approximations based only on contact stress and rolling diameter (ref. 8). Be- 
cause of the similarity in the failure mechanism between traction drive components and 
those experienced with rolling-element bearings, it is anticipated that the bearing fatigue 
life theory of Lundberg and Palmgren (ref. 9) can be successfully adapted to predicting 
traction drive life , In a similar vein, the work reported in reference 10 applied the 
statistical approach used by Lundberg and Palmgren for rolling -element bearings to life 
predictions for spur and helical gears. The surface fatigue life predictions from that 
analysis showed good agreement with actual test data. 

In view of the aforementioned, the objective of the work reported herein is to (1) de- 
velop an analysis based on classical rolling- element bearing fatigue theory for predicting 
the design life of traction drive systems and (2) to demonstrate the use of this analysis 
by forecasting the fatigue life expectancy of an advanced toroidal traction transmission 
vinder a variety of operating conditions. 


SYMBOLS 


a semimajor axis, m (in.) 

b semiminor axis, m (in.) 

c orthogonal shear stress exponent 

E Young's modulus, N/m (psi) 

S elliptic integral of the first kind 

e Weibull exponent 

F curvature difference 

axial force, N (lb) 

^ elliptic integral of the second kind 
g defined by eq. (18) 

H life, hr 

h depth to critical stress exponent 
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K,Kj constants of proportionality 

k ellipticity ratio 

L life (millions of input shaft revolutions) 

I length of rolling track, m (in.) 

m drive ratio, (output speed/input speed) 

n speed, rpm 

actual roll body load, N (lb) 

Qjj required roll body load, N (lb) 

0 

q maximum Hertz stress, N/m (psi) 

R cam radius, m (in.) 

r radius, m (in.) 

input torque, N-m (in. -lb) 
u stress cycles per revolution 

O O 

V stressed volume, m (in. ) 

w semiwidth of rolling track, m (in.) 

number of cavities 
number of rollers per cavity 

Zq depth to maximum orthogonal reversing shear stress 

a roller tilt angle , rad 

preload cam angle, rad 

0 roller cone angle, rad 

p traction coefficient 

1 Poisson’s ratio 

p curvature sum, m“^ (in. "^) 

Tq maximum orthogonal reversing shear stress, N/m^ (psi) 

i/' overload ratio 

Subscripts: 

A,B elastic bodies 

act actual 

c cavity 
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d drive 

in input contact 

o output contact 

r roller 

req required 
t torus 

x,y reference planes 


ANALYSIS 

Fatigue Life Analysis 

In 1947 Lundberg and Palmgren (ref. 9) published a theory for the failure distribu- 
tion of ball and roller bearings. The mode of failure was assumed to be subsurface 
originated (SSO) fatigue pitting . Lundberg and Palmgren theorized that SSO fatigue 
pitting was due to high stresses in the neighborhood of a stress raising incongruity in 
the bearing material. The important parameters are number of stress cycles L, mag- 
nitude of the critical stress amount of volume stressed V, and depth below the sur- 
face at which the critical stress occurs z^. The theory is widely used to predict 
rolling -element bearing fatigue life and was recently adapted for predicting the life of 
spur and helical gears (refs. 11 and 12). In reference 11 the formulation for life pre- 
diction of a steel gear set was confirmed with life data from full-scale spur gear tests. 
For a steel rolling element the number of stress cycles endured before failure occurs is 
given by the following equation (ref. 12); 


L = 



( 1 ) 


This equation is a modified form of the Lundberg -Palmgren theory for contact-fatigue 
life prediction and is applicable to gears, bearings, and other rolling- contact elements. 
For rolling -element bearings (and bodies in rolling contact in general) made of AKI 
52100 steel, Rockwell C62 hardness, with bearing life at a 90 percent probability of s\ir- 
vival, the following values are appropriate for use in equation (1) to determine life in 
millions of stress cycles (ref. 12): 


5 


I 



( 2 ) 


Kj = 1.428x10®^ (N and m units) 

= 3. 583x10^® (lb and in. units) 

e = — (for elliptical shape point contact) 

9 

= — (for line contact) (3) 

2 

h = 2| (4) 

c = lOj (5) 

The stressed volume is given by 

V = (6) 

where I is the length of the rolling track which is traversed during one revolution. The 
semiwidth of the rolling track is designated w, and is the depth to the maximum 
orthogonal reversing shear stress. 

If the body in question is subjected to several different contact loads during the load 

cycle, then Miner’s rtale (ref. 13) is useful in determining the element life as follows. 

Assume that loads Qj (j = 1 to k) act during one load cycle of body i. First, each of 

the lives L.. (j = 1 to k) is calculated by using equation (1) and assuming only the ap- 
13 th. 

propriate load j to be acting. Then the life of the i body when subjected to all k 

loads during a cycle is given by 



Heretofore the lives of the various elements were given in terms of millions of 
stress cycles. All bodies in the drive accumulate stress cycles at different rates be- 
cause their speeds of rotation and number of stress cycles per revolution are not all the 
same. In order to compare lives of the various bodies clock time should be used. As- 
sume that the speed in revolutions per minute of the i''*^ body is n^ and that there are 
Uj stress cycles per revolution. Then the life of body i in hours is given by 
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( 8 ) 


The life of the drive system is then found by applyii^ Weibull’s rvile (ref. 15) . If the 
drive system consists of n roll bodies and the life of each is designated (i = 1 
to n) , then the drive life is given by 


«d = 



(9) 


Contact Stress Analysis 


The stress analysis of elastic bodies in contact was developed by H. Hertz (ref. 14) . 
Hertz assumed homogeneous solid elastic bodies made of isotropic material which are 
characterized by Young’s modulus E and Poisson’s ratio Bodies A and B in con- 
tact are assumed to have quadratic surfaces in the neighborhood of the contact point. 

Figure 1 shows two bodies in contact. Planes x and y are the respective planes 
of maximum and minimum relative cxmvature for the bodies. These planes, called the 
principal planes, are mutuaUy perpendicular. Planes x and y must be chosen so that 
the relative curvature in plane x is greater than in plane y : 



^Ax ^Bx ^Ay ^By 


( 10 ) 


The radii of curvature may be positive or negative depending on whether the surfaces 
are convex or concave, respectively. 

After the bodies are pressed together the contact point is assmned to flatten into a 
small area of contact which is bounded by an ellipse with major axis 2a and minor 
axis 2b as shown in figure 1 . Plane y contains the major axis of the contact ellipse 
and plane x contains the minor axis. The ratio k = a/b is called the ellipticity ratio 
of the contact. The values of k range from 1 to for various curvature combina- 
tions of contacting surfaces. For cylinders in contact, the ellipticity ratio is <», and 
the flattened area of contact is a rectangular strip. For spheres in contact the eUiptic- 
ity ratio is 1 . The first type is called line contact and aU other types are called point 
contact. The theory of Hertz is summarized by Harris (ref. 15). The relation between 
ellipticity ratio and the geometry of the contactii^ bodies is given by the following 
transcendental equation: 
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( 11 ) 


F = 



P 


(k^ + 1) -2^ 

(k^ - l)i 


p 1 , 1 , 1 , 1 

^Ax ^Bx ^Ay ^By 


/ n/2 


|1 _ (1 - J-jsin^<p 

k2 


■ 1/2 


d(p 


I 


v/2 


1 - [l — i-]sin2^ 


1/2 




( 12 ) 


(13) 


(14) 


where F is the ciorvature difference, p the curvature sum, and S and ^ are elliptic 
integrals of the first and second kind, respectively. Hamrock has developed an iterative 
scheme by which these equations may be solved (ref. 16). Hamrock’s method was used 
to perform the stress calculations needed for this analysis . 

The maximum surface contact pressure at the center of the ellipse is 


q = 


3Q 

2ffab 


where 


a = a*g 


b = b*g 



(15) 


(16) 

(17) 


(18) 


(19) 
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( 20 ) 


b* = 



The maximum reversing orthogonal shear stress occurs at a depth imder- 
neath the surface of contact: 
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V2t - 1 

2t(t + 1) 


q 


1 

(t + l)\/2t - 1 


b 


where t is an auxiliary parameter which is related to k by 


k = 


[(t^ - l)(2t - 1) 


- 1/2 


( 21 ) 

( 22 ) 

(23) 


Table I may be used to calciolate the stresses and contact ellipse dimensions when it is 
impractical to solve equation (11) . 


Analysis of Toroidal Drive 

A schematic diagram of the single cavity cone-roller toroidal (CRT) drive is shown 
in figure 2. This drive is of the same contact geometry as that presented in reference 7. 
The drive consists of one input and one output traction disk. The drive disk conforms 
to a toroidal shaped cavity with radius r^. The radius of the torus is r^. The output 
disk is driven by an input disk through the intermediate rollers, which are crowned with 
radius r^. The intermediate roller is compressed between the drive disks by the pre- 
load cams in order to develop the contact loads necessary for transmitting the required 
power without gross slippage. 

The cone rollers are offset from the cavity center by the amount r^ cos 0/2. The 
’ ai^le of tilt a determines the drive ratio. For positive values of a tiiere is reduction 
in speed, and for negative values of a the speed is increased through the drive. The 
• equation relating drive ratio m to tilt angle a is 
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m = 


(24) 


^in 


r 


o 


^t" 


^t- 



Figure 3 shows a schematic of the dual cavity CRT drive. This dual cavity version 
is identical in contact geometry to the single cavity drive . Two cavities eliminate the 
need for large input and output thrust bearings. Notwithstanding the added size and 
complexity, it is probable that the dual cavity arrangement will be the first choice for 
an automotive application. 

The traction drive is able to transmit power due to the resistance to slip that is 
developed in the thin film of oil in the rolling contact. The tractive force is a function 
of the type of oil, temperature, rolling and sliding speeds, degree of spin, surface 
finish, and normal contact force per unit width of contact. In order to prevent slipping 
at the traction contacts the normal load must be maintained above a certain minimum 
level so that the tractive force provided by the oil exceeds that required by the contact. 
This normal load is maintained by a rolling cam preload mechanism as shown in fig- 
ure 4. The axial thrust F„ is related to the input torque T. by 


^a = 


(25) 


where 


K = 


1 


R tan /3 


(26) 


The ayial force then acts to pinch the cone roller and thus apply the necessary normal 
force. The actual normal force produced is given by 


Q 




act 


sini 


(i-a 


V2 


x„x. 


c r 


sin 


/e 

- - Q! 

\2 


^c^r 


(27) 


This normal force produced is only a function of the drive ratio, cam constant, and 
cone angle. However, the normal force that is required to assure sufficient traction is 
a function of all the factors previously mentioned that cause the formation of the EHD 
film at the disk -roller contact points. The traction coefficients can range from 0.03 to 
0.08 (ref. 7). However, under the operating conditions which exist in the CRT, the 
practical upper limit is about 0.06. In the present analysis a constant available traction 
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coefficient of 0.06 is assumed at all conditions. This assumption allows the most gen- 
erotis relation between traction and compressive normal forces . Therefore the life 
analysis is optimistic for this type of drive. 

The overload factor is the ratio of the actual normal force on the contact to the 
normal force required to transmit the given torque. The required contact force is 
given by 


Q. 


T. 

in 


req 


fir 


in^c^^r 


(28) 


where 


From equations (27), (28), 


^in ^t ■ ^c 


(j'“) 


and (29) the expression for overload factor is obtained; 




^act fi 

Q tan B 
^req ^ 




(29) 


(30) 


The minimum overload occurs at roller tilt angle given by 

(31) 

To provide some margin of safety against gross slippii^, the minimum overload factor 
should be significantly larger than unity. In the present analysis a minimum overload 
factor of 1 . 2 was assumed . 

The contact radii that are needed for the stress analysis may be found from figure 2. 
At the constant of input disk and roller the principal radii of curvature are 



T 



^B1 = ^c 


(32) 


(33) 
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1 



(34) 


^A2 = -^c 


^B2 = 


For the contact of output disk and roller the radii are 






COS 


^Bl=^c 


^A2 = -^c 


^B2 ■ ^r 


(35) 

(36) 

(37) 

(38) 

(39) 


In order to determine the orientation of the contact ellipse in relation to the rolling track 
the inequality of equation (10) must be satisfied. Since it is not known a priori whether 
plane 1, the plane of rolling, or plane 2, the transverse plane, is the plane of the max- 
imum curvature, the subscripts x and y were not used in the previous general eqvia- 
tions. 

If reference plane 1 is the plane of maximiun curvature (plane x) according to 
equation (1), then the major axis of the contact ellipse is perpendircular to the rolling 
track. Conversely, the minor axis of the contact ellipse is perpendicular to the rolling 
track if plane 2, the transverse plane, is the plane of maximum curvature (plane x). 

For the dimensions and drive ratios used herein the plane of maximum curvature is in- 
deed plane 1 and hence the major axis of the contact ellipse defines the width of the 
rolling track. This is true for both the input and output traction contacts. 

The 90-percent life for each component of the traction drive is found for a given 
drive ratio by applying equations (1) and (7) . The element lives given by equations (1) 
and (7) are in terms of stress cycles endured. The lives may also be expressed in hours 
of operation before pitting damage occurs by using equation (8) . The input disk endures 

stress cycles per revolution if there are roller elements. The life of the input 
disk is given by 


H = 

'”'Vin\ 60 / 


(40) 


I 


♦ 


The life of the output disk is given by 
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(41) 


H 


o 


r o 






The intermediate roller elements contact both the input and output disks. For a unity 
drive ratio, the stress and contact ellipse dimensions are the same at both contacts. In 
this case the roller is subjected to two equal stress cycles per roller revolution. 

When the drive ratio is not unity, the roller is subjected to two different states of 
contact stress per revolution. In this case, the life of the roller is denoted by (L^Hn 
when the input is assumed to act alone on the roller. Similarly, (Lj.)q is the life if 
only the output contact load is acting. Then using Miner's rule (eq. (7)) the expected 
life of the roller is calculated: 

-1 

(42) 


Lr = 




(LPo 


In hours, the life of a single roller is determined by equation (8) : 


^ n^[ 60 j n.„ r.„ Uo 


(43) 


The life of the single cavity drive consisting of the manifold of four elements (2 drive 
disks plus 2 rollers) is found according to Weibull's rule (eq. (9)): 


H 


d 




(44) 


RESULTS AND DISCUSSION 

A single cavity CRT drive was analyzed for life and Hertz contact stress. The di- 
mensions of the drive were taken eqvial to those of the drive in reference 7, and are as 
follows : 
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Torus diameter, cm (in.) 11.43 (4.50) 

Cavity diameter, cm (in.) 8.84 (3.48) 

Roller crown radius, cm (in.) 3.12 (1.23) 


Figure 5 gives the overload factor for the drive as a function of ratio from equa- 
tion (29). The minimvim overload factor occurs at drive ratio 0.647. This minimum 
overload factor was set equal to 1 . 2 to provide a 20 percent margin of safety against 
slip. The variation in the overload factor with ratio is primarily dependent on traction 
drive geometry. If the overload factor becomes too large then the fatigue life of the 
drive will be penalized. The cone roller toroidal drive under study experiences signif- 
icantly less variation in the overload factor than other comparable toroidal roller drives 
(ref. 7). For the case presented, the overload factor ranged from 1.2 to 1.8. 


Parametric Variations 

Figure 6 summarizes the results of the stress analysis. The analysis showed the 
Hertz stress at each traction zone is proportional to the cvibe root of input torque for a 
constant drive ratio. If the torque is held constant, there is a 2. 3 to 1 variation in the 
Hertz stress at the input disk -roller contact as the drive ratio increases from 1:3 to 
4:1. The larger maximum Hertz stress is on the input contact for drive ratios smaller 
than 1:1 and on the output contact for higher ratios. 

Figure 7 shows the load-life curves that resulted from the analysis. Three differ- 
ent ratios were considered. The curves presented show that for any fixed operating ra- 
tio the life is related to load by 


Lcc(TJ-> 

The most severe service condition for any given load is the 1:3 reduction ratio. 

The lives of the individual components were also calculated. The relative lives of 
the drive disks and intermediate roller are given in table II. The lives are normalized 
by the drive system life for each ratio considered. For low ratio the life limiting com- 
ponent is the input disk and for the high ratio the life is limited by the output disk. If 
the component lives are normalized by the system life at 1:3 drive ratio, then the effect 
of ratio on each component life is included. Those results are in table HI. 

Table IV is a summary of the effect that geometry changes have on life and Hertz 
stress. Life and stress calculations were performed for scaled up versions of the basic 
single cavity CRT drive with two rollers. Also, the effect of adding more rollers, as 
well as one more cavity, was investigated. It was found that by a size increase of 10 per- 
cent a life equivalent to the three-roller drive (case 5) was obtained. To achieve life par- 
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ity with the double cavity drive the basic drive size had to be increased by about 20 per 
cent. The drive life increases rapidly with size according to the proportionality 

O A 

L oc (size factor) ’ 

Therefore, if the costs of manufacture are considered it is more cost effective to sim- 
ply increase size in order to achieve higher load ratings. However, with increased 
sizes the rolling speed at the traction contacts increases and the contact stress de- 
creases. It is known that the traction coefficient decreases if rolling speed increases 
or contact stress decreases (ref. 17). Therefore, due consideration must be given to 
providing sufficient preload to compensate for the decreased traction coefficient that is 
available. A derating of the life expectation will be the final result of such a compensa- 
tion adjustment. 


Design Study 

A design study similar to that presented in reference 7 was also performed. The 
double cavity drive (case 6) was incorporated into the vehicle drive train as shown in 
figure 8. This drive system is modeled after the system described in the aforemen- 
tioned reference which used the 4 -liter (250 CID) six-cylinder engine and Plymouth B 
body. For the range of driving conditions considered, the torque input to the toroidal 
drive ranged up to 121 newton-meters (1071 in. -lb), which corresponds to an engine out- 
put of 102 horsepower at 4000 rpm. Figure 9 shows the calculated life as a function of 
drive ratio. Figure 10 shows the calculated values of maximum Hertz stress as a func- 
tion of drive ratio. 

For combinations of low ratio and high torque due to sudden overloads, the maxi- 
mum Hertz stress on the input contact can exceed the AFBMA brinell limit of 517 000 
newtons per square centimeter (750 000 psi). The brinell limit is defined as that maxi- 
mum Hertz stress which would produce identations of 0.0001 centimeter per centimeter 
of roller diameter of AISI 52100 steel of Rockwell C64 hardness (ref. 18). In refer- 
ence 19, figure 7 also shows that a Hertz stress of 517 000 newtons per square centime- 
ter (750 000 psi) will cause an indentation of 0.0001 times the ball diameter. However, 
the elastic limit of the material is exceeded at much lower stresses, and the 0.0001 
times diameter criterion may not apply as a definition of excessive deformation for a 
toroidal drive . This is because the rolling elements in bearings generally run over the 
same track; therefore, a slight amount of grooving would have little adverse effect. How- 
ever, since the intermediate rollers in a toroidal drive are required to sweep over the 
surface of the input and output disks in order to change ratio, even a small amount of 
grooving can have a detrimental effect on roller component life and the ability to adjust 


15 



drive ratio. In reference 15, equations for plastic deformation of ball bearings are given 
in chapter 11, Using those equations, stresses of only 345 000 to 414 000 newtons per 
square centimeter (500 000 to 600 000 psi) are required to cause a deformation of 0. 0001 
times the ball diameter. 

In order to make the results of the analysis more useful for design assessment, cal- 
culated values of life and maximum Hertz stress are plotted against vehicle speed in fig- 
ure 11. The life is given in hours for three different engine conditions of power output 
and speed. Also shown as dashed lines are the tire slip limit and the zero grade - zero 
wind road load conditions which span the probable values of vehicle operating parameters. 
If a representative vehicle operating point of 48 kilometers per hour (30 mph), 29 engine 
horsepower at 2000 rpm, and 1. 2 drive ratio is selected, the resulting theoretical life is 
331 hours at a 90-percent probability of survival. This condition can be taken as an av- 
erage approximation of the operating parameters over the service life of the vehicle. If 
the vehicle is operated continuously at the tire slip limit, the drive life would be just un- 
der 10 hours. 


Life Adjustment Factors 

In recent years, improvements in rolling-element bearing analysis, materials, 
processing, and manufacturing techniques have generally permitted an increase in the 
expected life of a bearing contact for a given application. In reference 20 bearing 
life adjustment factors have been developed which, when applied to Lundberg-Palmgren 
fatigue life predictions, will provide much better estimates of bearing fatigue perform- 
ance. Many of these factors are also applicable to toroidal drive life predictions. The 
factors considered are (1) materials factor, (2) processing factor, and (3) lubrication 
factor . 

It is anticipated that the contacting elements of the CRT drive will be made from 
consumable vacuum melted (CVM) AISI 52100 steel which has become the bearing in- 
dustry's standard material. According to reference 20, the use of CVM processed AISI 
52100 steel in a bearing application would result in a life amplification factor of 6 to be 
applied to theoretical life estimates . 

With regard to traction lubricant fatigue life effects, some investigators have sug- 
gested that certain types of traction lubricants (viz, synthetic cycloaliphatic fluids) may 
also show superior fatigue life performance when compared to conventional lubricants 
under identical operating conditions (ref. 6) . Unfortunately, the life improvement po- 
tential of these types of traction oils is somewhat in dispute in view of the research of 
reference 21, which could not establish any statistically significant fatigue life advan- 
tages. At this time there is no firm basis for assigning life improvement factors based 
on the traction lubricants' chemical effects. 
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However, it is well known that the thickness of lubricant film which separates con- 
tacting machine elements can have a strong influence on the contacts’ fatigue life. In 
reference 20, the effectiveness of a lubricant film in terms of the ratio of film thickness 
to composite surface roughness h/a is related to fatigue life. The lubricant -fatigue 
factor has been assigned a value of 1.0 when h/a has a value of approximately 1.3. 

In the case of the CRT drive for the representative operating point of 29 horsepower, 
48.3 kilometers per hour (30 mph), and a drive ratio of 1.2, the calculated film thick- 
nesses of the input and output contacts are 0. 51 and 0. 48 micrometer (20 and 19 pin.), 
respectively, using the Archard and Cowking film formula (ref. 20). This calculation, 
the details of which appear in the appendix utilizes an advanced commercially available 
traction fluid with a contact inlet temperature of 344 K (160° F) . These film thickness 
values must be regarded as quite marginal considering that they were calculated for a 
rather modest operating condition and that a surface finish of 0.61 micrometer (24 pin.) 
rms has been proposed for the CRT drive elements (ref. 7). In fact, the calculated in- 
put and output h/a values of 0.59 and 0. 56, respectively, are so low at this condition 
that it is questionable whether long-term operation can be sustained without drastic sur- 
face distress. If operation is assumed possible, an 80-percent reduction in fatigue life 
would be expected (ref. 20). Thus, the sixfold fatigue life increase due to the CVM 
52100 steel is offset by a lubrication factor of 0.2 for a net 1.2 life improvement factor. 

In summary, the expected life at the selected duty cycle point based on classical 
bearing fatigue theory and the experience factors mentioned is 397 hours (1.2 x 331). 

The expected life reported in reference 7 for the identical conditions is approxi- 
mately 38 times longer. The mileage life using the methods reported herein is approxi- 
mately 19 200 kilometers (11 900 miles) for this representative condition. Based on the 
life-scaling relations presented earlier, the CRT drive size should be scaled up to ap- 
proximately 130 percent of its current size to provide for 160 900 kilometers (100 000 
miles) of service. This approach is acceptable if the anticipated loss in the available 
traction coefficient due to increased surface speeds and lower contact pressures is 
tolerable . 

It should be pointed out that at operating conditions which are more stringent than the 
condition examined, insufficient lubricant film will be generated to separate the relatively 
rough 0.61 micron (24 pin.) rms roller svirfaces and gross surface damage wovild most 
likely result. To combat loss of lubricant film the speed of the contact should be in- 
creased, the temperature of the lubricant shoxild be lowered, and the surface finish of 
the roller elements should be greatly improved. These improvements which will un- 
doubtedly add to the cost and possibly the size of the drive are necessary for reliable 
operation. 
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SUMMARY OF RESULTS 


A contact fatigue life analysis for traction drives has been developed. The life 
analysis was based on a modification of Lundberg-Palmgren theory which is the basis 
of rating bearing life in the rolling-element bearing industry. 

The fatigue life analysis was used to do a parametric study of a toroidal type trac - 
tion drive . The design study was compared to previously published work for the dual 
cavity toroidal drive as applied to a typical compact passenger vehicle drive train. 

Drive life and contact stresses were computed using a computer program based on the 
method presented herein. A constant traction coefficient of 0.06 was selected and a 
minimum contact overload factor of 1.2 to preclude slip was used. Conditions of elasto- 
hydrodynamic lubrication of the contacts were computed at a selected operating point and 
life adjustment factors due to advances in materials and lubricants technology were con- 
sidered. 

The following results were obtained from a study of a toroidal drive with a torus di- 
ameter of 11. 43 centimeters (4. 50 in. ), a cavity diameter of 8. 84 centimeters (3. 48 in. ), 
and a roller crown radius of 3. 12 centimeters (1, 23 in. ); 

1. For a representative duty cycle condition wherein the engine delivers 29 horse- 
power at 2000 rpm with the vehicle moving at 48.3 kilometers per hour (30 mph) at a 
drive ratio of 1.2, the life of the CRT drive is 19 200 kilometers (11 900 mile). This is 
38 times less than previously stated in the literature for the same conditions . 

2. For the previously stated conditions and a surface finish of 0.61 micrometer 
(24 pin.) rms, the ratio of lubricant film thickness to composite siirface roughness was 
less than 0.6. The lubricant film is barely sufficient to prevent severe roller component 
surface distress. 

3. For conditions of high reduction ratio and high transmitted loads, the brinelling 
limit of high quality, hardened bearing steel is approach for the input disk -roller con- 
tact. Shock or dynamic overloading may cause significant brinelling of the input disk or 
roller . 

4. Life is directly proportional to the 8.4 power of size. To provide for 160 900 
kilometers (100 000 mile) of vehicle service life at the aforementioned duty point the CRT 
drive size must be scaled to approximately 130 percent of its ciirrent size. This can 

be accomplished only if the loss in available traction coefficient due to the increased sur- 
face speeds and reduced contact pressures is not too severe. 
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5. Drive life is inversely proportional to the cube of torque. 

6. Hertz stress is proportional to cube root of input torque. 

Lewis Research Center, 

National Aeronautics and Space Administration, 
and 

U.S. Army Air Mobility R&D Laboratory, 

Cleveland, Ohio, July 9, 1976, 

505-04. 



APPENDIX - ELASTOHYDRODYNAMIC FILM THICKNESS CALCULATION 


The following lubricant film thickness equation for elliptical contacts was developed 
by Archard and Cowking (as reported in ref. 20): 
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where 

h 


Ui,U2 

Q 


E 

P-Q 

a 


film thickness, m (in.) 

surface velocities of rolling elements, m/sec (in. /sec) 
contact load between rolling elements, N (lb) 

equivalent contact radii in orthogonal directions, 1 refers to direction of 
rolling, m (in.) 

2 

Young's modulus of elasticity, N/m (psi) 

Poisson's ratio 

2 2 

lubricant absolute viscosity at atmospheric pressure, N-sec/m (Ib-sec/in. ) 

2 -1 

reciprocal asymptotic isoviscous pressure -viscosity coefficient, (N/m )" , 
(psi‘^) 


The equivalent radius of contacting bodies in the plane of rotation is given by 
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Ri 
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and in the plane transverse to the direction of rolling by 

.-1 


^2 ~ (^Ag 


O' 


(A3) 


where ’ ^B ’ ^B principal radii of curvature as defined before. 

12 1 2 

At the representative duty cycle point the following conditions exist: 
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Vehicle speed, km/hr (mph) 48.78 (30) 

Engine power, hp 29 

Engine speed, rpm 2000 

Drive input speed, rpm 3000 

Drive ratio 1.1653 

Contact load between rolling elements, N (lb) 11 782 (2648.8) 

Surface velocities of rolling elements, Uj^ = U 2 , m/sec (in. /sec) 9.739 (383.43) 

Young’s modulus of elasticity, E, N/m^ (psi) 2.07x10^^ (30x10®) 

Poisson's ratio, | 0.3 

Equivalent contact radii in orthogonal directions, m (in.): 

(input) 2. 397 x10"2 (0.9438) 

Rg (input) 10 . 66X10'^ (4 . 1965) 

Rj (output) 2. 057X10"^ (0.8099) 

R 2 (output) 10.66x10“^ (4. 1965) 

For a synthetic cycloaliphatic traction fluid at a temperature of 344 K (160° F), 

|i =0.010 N-sec/m^ (1.452x10'® (Ib-sec/in. ^)) (ref. 21) (A4) 


a = 2.06x10'® m^/N (1.42x10"'^ psi'^) (ref. 22) (A5) 


The following film thicknesses were calculated for the input and output contacts: 


‘'input = (20x10'® in.) 

(A6) 

'output 

(A7) 


Usually the root mean square (rms) surface finishes of the contacting bodies 
and Og are used to determine the composite surface roughness as follows: 


{ 2 

a = ^Oi + 


(A8) 


In reference 7 the proposed roller surface finishes are 0.61 micrometer (24 pin. 
rms); thus, ct = 0.86 micrometer (33.9 pin.) and the h/a ratios for the input and output 
are 



(A9) 


(AlO) 
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TABLE I, - VALUES OF CONTACT 


STRESS PARAMETERS 


F 

a* 

b* 

‘ 1 

0 

1 

1 

1.2808 

.1075 

1.0760 

.9318 

1.2302 

.3204 

1.2623 

.8114 

1.1483 

.4795 

1.4556 

.7278 

1.0993 

.5916 

1.6440 

.6687 

1.0701 

.6716 

1.8258 

.6245 

1.0517 

.7332 

2.011 

.5881 

1.0389 

.7948 

2.265 

,5480 

1.0274 

.83495 

2.494 

.5186 

1.0206 

. 87366 

2.800 

.4863 

1.0146 

.90999 

3.233 

,4499 

1.00946 

.93657 

3.738 

.4166 

1.00612 

.95738 

4.395 

.3830 

1.00376 

.97290 

5.267 

.3490 

1.00218 

.983797 

6.448 

' .3150 

1.00119 

.990902 

8.062 

.2814 

1.000608 

.995112 

10.222 

.2497 

1.000298 

.997300 

12.789 

.2232 

1.000152 

.9981847 

14.839 

.2070 

1.000097 

.9989156 

17.974 

.18822 

1.000055 

.9994785 

23.55 

. 16442 

1.000024 

.9998527 

37.38 

.13050 

1 . 000006 

1 

OO 

0. 

1.000000 


TABLE II. - COMPONENT LIVES NORMALIZED BY 
DRIVE SYSTEMS LIFE AT VARIOUS RATIOS 


AND CONSTANT LOAD 


Ratio 

— 

Drive 

system 

Input disk 

Output disk 

Single roller 

1:3 

1 

1.75 

371 

3.75 

1:1 

1 

3.45 

3.45 

3,51 

4:1 

1 

1498 

1.73 

3,78 
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TABLE ni. - COMPONENT LIVES NORMALIZED BY 


DRIVE SYSTEM LIFE AT CONSTANT 


LOAD AND 1:3 RATIO 


Ratio 

Drive 

system 

Input disk 

Output disk 

Single roller 

1:3 

1 

1.75 

371 

3.75 

1:1 

43.6 

151 

151 

153 

4:1 

2.84 

4253 

4.92 

10.73 


TABLE IV. - EFFECT OF SIZE, NUMBER OF ROLLERS, AND 
NUMBER OF CAVITY BANKS ON CONTACT FATIGUE 


LIFE AND HERTZ STRESS 


Drive configuration 

For any fixed ratio and load 


Life increase, 
percent 

Reduction in 
Hertz stress, 
percent 

1. Basic CRT drive: 

11.43 cm (4.50 in.) torus diameter 



8.84 cm (3.48 in.) cavity diameter 
3. 12 cm (1.23 in.) roller crown ra- 
dius single cavity with 2 rollers 
2 . 10 Percent larger 

123 

9 

3. 20 Percent larger 

363 

17 

4. Double size 

34 000 

50 

5. Same as 1 except 3 rollers 

129 

13 

6. Same as 1 except 2 cavities with 

329 

21 

2 rollers per cavity 







rEdge of plane 1 (plane of rolling) 



Figure 2. - Single cavity cone- roller toroidal drive. Ratio continuously variable by tilting 
plane of rollers. Tilt designated by angle a. Drive rollers shown in speed reducer posi- 
tion. 


26 







Input shaft torque, in. -lb 
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